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Abstract — We are concerned with catalyst-assisted probabilis- 
tic entanglement transformations. A necessary and sufficient 
condition is presented under which there exist partial catalysts 
that can increase the maximal transforming probability of a given 
entanglement transformation. We also design an algorithm which 
leads to an efficient method for finding the most economical 
partial catalysts with minimal dimension. The mathematical 
structure of catalyst-assisted probabilistic transformation is care- 
fully investigated. 

Index Terms — Probabilistic entanglement transformation, ma- 
jorization, catalysis, catalyst-assisted transformation, partial cat- 
alyst. 



I. Introduction 

QUANTUM entanglement plays an essential role in quan- 
tum information processing [1]. Indeed, it is a necessary 
resource in quantum cryptography [2], quantum superdense 
coding [3], and quantum teleportation [4], which are striking 
tasks in quantum information processing. When entanglement 
is treated as a type of resource, the study of how to quantify 
and manipulate it becomes crucial (for a survey of quantum 
information theory, we refer to [5]). A fruitful research di- 
rection is to try to discover the laws that must be obeyed 
when transforming between different forms of entanglement 
using only local operations on the separate subsystems and 
classical communication between them. This kind of transfor- 
mation is usually abbreviated as LOCC. The communication 
constraints that characterize LOCC are fundamentally and 
practically important, since many applications of quantum 
information processing involve spatially separated parties who 
must manipulate an entanglement state without performing 
joint operations. 

Suppose two spatially separated parties, Alice and Bob, 
share a bipartite quantum state \tpi) £ C" (X)C" with Schmidt 
decomposition 



want to transform this initial state into a desired bipartite state 
\-ip2) with Schmidt decomposition 



IV'i 
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where ai > . . . > a„ > are the Schmidt coefficients of 
IV'i) and ^^ai = 1. and jis) are orthonormal bases of 
Alice's and Bob's systems, respectively. Suppose the parties 
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where /3i > . . . > /3n > 0, and X)i A = 1- The orthonormal 
bases and (also and lis)) are not necessarily 
the same. Nielsen [6] proved that Alice and Bob can realize 
this transformation from \'4>i) to |-02) by LOCC if and only if 



I 

ai < ^ (3i for any 1 < ^ < n, 



i=l 



i=i 



with equality holding when I = n, ot equivalently, by the 
theory of majorization [7] [8], A^j is majorized by A^j, written 

A^i -< A^2 1 

where the probability vectors A^^ and X^^ denote the Schmidt 
coefficient vectors of \tpi) and |'02), respectively. 

Nielsen's work establishes a connection between the theory 
of majorization in linear algebra and entanglement transfor- 
mation. Furthermore, since the necessary and sufficient condi- 
tion mentioned above is very easy to check, it is extremely 
useful in telling whether one bipartite entangled state can 
be transformed into another by LOCC. Nielsen's theorem 
directly implies that there exist incomparable states in the 
sense that any one cannot be transformed into another only 
using LOCC. To treat the case of transformations between 
incomparable states, Vidal [9] generalized Nielsen's work by 
allowing probabilistic transformations. He found that although 
deterministic transformation cannot be realized between in- 
comparable bipartite pure states, a probabilistic one is always 
possible (notice that when multipartite states are considered, 
this statement does not hold [10]). Furthermore, he gave an 
explicit expression of the maximal probability of transforming 
one state to another To be more specific, let P{\ipi) IV'2)) 
denote the maximal probability of transforming into \ip2) 
by LOCC. Then 
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where Ei{\^^) denotes ^i- what follows, we extend 
this notation to any probability vector 

Another interesting phenomenon discovered by Jonathan 
and Plenio [11] is that sometimes an entangled state can enable 
otherwise impossible entanglement transformations without 
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being consumed at all. A simple but well known example is 

^ |f/'2) but IV'i) ® \4>) IV'2) ® 10), where 

1^1) = VO|00) + VqA\11) + \/0l|22) + \/0l|33), 
1^2) = \/a5|00) + \/a25|ll) + \/a25|22), 

and 

=ya6|44) + Va4|55). 

The role of the state |(^) in this transformation is analogous 
to that of a catalyst in a chemical process. The mathe- 
matical structure of this phenomenon, so called catalyst- 
assisted entanglement transformation, was carefully examined 
by Daftuar and Klimesh [12]. They found that there does 
not exist an upper bound on the dimension of catalysts that 
should be considered, in trying to determine which states can 
be transformed into a given state. Furthermore, they proved 
that any nonuniform state, which has at least two nonzero 
Schmidt coefficients nonequal, can serve as a catalyst for some 
entanglement transformation. On the other hand, Eisert and 
Wilkens found that catalysis is also helpful in entanglement 
transformations for bipartite mixed states [13]. 

In this paper, we examine the power of catalysis in proba- 
bilistic entanglement transformations. We have noticed that 
in [11] Jonathan and Plenio found that in some cases, an 
appropriately chosen catalyst can increase the maximal trans- 
formation probabiUty of incomparable states. The example 
they presented is as follows. Let 

IV^i) = \/a6|oo) -I- \/a2|ii) + \/a2|22) 

and 

= V(15\00) + VoA\n) + VoJ\22). 

The maximal probability of transforming into IV'2) under 
LOCC is 80% while when a catalyst 

10) = V0ii5|33) + \/(I35|55) 

is introduced, the probability can be increased to 90.4%. 
They also showed that enhancement of the transformation 
probability is not always possible. However, no further results 
on this topic were given in their paper. 

The main aim of our paper is to study the structure of catal- 
ysis as applied to probabilistic entanglement transformations. 
We give a necessary and sufficient condition for the existence 
of partial catalysts (quantum states which can increase the 
maximal transforming probability while not being consumed) 
for a given entanglement transformation. Rather surprisingly, 
we find that whether or not the probabiUty can be increased 
depends only on the minimal Schmidt coefficients of the 
original state and the target state, provided that the maximal 
transforming probability is less than 1 (the probabiUty cannot 
of course be increased if equal to 1). To be specific, a 
probabilistic transformation has partial catalysts if and only if 
the maximal transforming probabiUty is less than the minimum 
of 1 and the ratio of the minimal Schmidt coefficient of the 
original state to that of the target state. Furthermore, we show 
that if the maximal probabiUty of a transformation can be 
increased by some catalyst, then there is a sequence of 2 x 2 



dimensional states that increases the maximal probabiUty of 
the transformation. 

For any given entanglement transformation, we present a 
systematic way to construct partial catalysts. The catalysts are, 
however, not economical in general in the sense that they 
do not necessarily have the minimal dimension among all 
partial catalysts for this transformation. In fact, the problem 
of constructing systematically the most economical partial 
catalysts for any given transformation seems to be hard and 
remains open. We can, however, give a numerical solution to 
this problem by solving a series of inequaUties. 

This paper is organized as follows. In Section II, we con- 
sider briefly probabilistic entanglement transformations with- 
out the aid of catalysis. We first provide a simple connection 
between probabiUstic transformations and deterministic ones, 
which is helpful in realizing probabilistic transformations since 
deterministic ones have been well researched. We then exam- 
ine properties such as monotonicity and continuity of the set of 
states that can be transformed into a given state by LOCC with 
a probability not less than a given positive number. Section III 
is the main body of this paper. We present here a necessary 
and sufficient condition under which a given probabilistic 
transformation has partial catalysts. Moreover, the catalysts 
are systematically constructed. To find the most economical 
ones, we present first an algorithm to decide whether there 
exist partial catalysts with a given dimension and find out 
all suitable ones. Based on this algorithm, the mission of 
constructing the most economical partial catalysts is achieved 
by applying the algorithm to state spaces with increasing 
dimension from 2 (an upper bound on the dimension we should 
consider can be predetermined because we have constructed a 
non-economical one). In Section IV, we generalize the result 
of Daftuar and KUmesh [12] to the set of states that can 
be transformed into a given state by catalyst-assisted LOCC, 
with a probability not less than a given positive number. We 
find that this set shares many properties with the well known 
set which consists of all states being trumped by the given 
state (for the latter set, we refer to [12] for details). To be 
more specific, the generaUzed set is convex and not closed in 
general; the dimensions of catalysts that should be considered 
in trying to determine the states in the set have no upper 
bound. We further investigate the mathematical structure of 
this generalized set and find out all the boundary and extreme 
points. This gives an answer to Nielsen's open problem in the 
case of deterministic transformations. FinaUy, a conclusion is 
drawn and some open problems are presented in Section V. 

For simplicity, in what follows we denote a bipartite quan- 
tum state by the probability vector of its Schmidt coefficients. 
This will not cause any confusion because it is well known that 
the fundamental properties of a bipartite state under LOCC are 
completely determined by its Schmidt coefficients. Therefore, 
from now on, we consider only probability vectors instead of 
bipartite quantum states and always identify a probability vec- 
tor with the corresponding quantum state. Sometimes we even 
omit the normaUzation of positive vectors to be probability 
ones for the sake of simpUcity. 
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II. Probabilistic Entanglement Transformation 

In this section, we discuss the structure of probabilistic 
entanglement transformations without catalysis. Denote by V" 
the set of all n-dimensional probability vectors and let x,y, . . . 
range over Given a positive number A < 1, let 

S\y) = {xeV'' ■.Pix^y)>X}, 

which is the set of n-dimensional probability vectors that can 
be transformed into y by LOCC with the maximal probability 
not less than A. When A = 1, the set reduces to the well 
known set S{y) which includes exactly the vectors that can 
be transformed into y with certainty, or equivalently, that are 
majorized by y (see [14] for details of S{y)). 

What we would like to point out first is that there is a 
simple relationship between probabilistic entanglement trans- 
formation and the theory of weak majorization, just like 
the connection between deterministic transformation and the 
theory of majorization discovered by Nielsen in [6]. Recall that 
an n-dimensional positive vector u is called super-majorized 
[7] by another n-dimensional positive vector v, written u 
V, if and only if 

n n 
i—l i—l 

for each / in the range 1 through n. Here denotes the 
vector obtained by arranging the components of u in nonin- 
creasing order Notice that the only difference between super- 
majorization and majorization is the omission of the equality 
requirement at ^ = 1. It is very easy to check by G. Vidal's 
formula that x e <S'^(?/), or equivalently, P{x ^ y) > A if 
and only if the super-majorization relation 

X Xy 

holds. 

It is well known that there is a close connection between 
majorization and doubly stochastic matrices [15]. To be spe- 
cific, for all x,y £ V^, x ^ y if and only if x = Dy for 
some doubly stochastic matrix D. Here a matrix D is called 
doubly stochastic if it is positivity preserving and every row 
and column sums to 1. That is, 

Vz, j : D,j > and Vj : ^ D,, =Y^Dj,^l. 

i i 

Unfortunately, to the authors' knowledge, super-majorization 
does not have such a correspondence. In order to make use 
of the known results about majorization, we must connect 
probabilistic entanglement transformation and majorization. 
The following lemma is just for this purpose. 

Lemma 1: For x,y E V"" and < \ < \, x E {y) if and 
only if X e S{y\), where 

yx = {l-XE-2{y),Xyi,---,\yL)- (D 

That is, S^iy) = Siyx). 

Proof: It follows directly from the definitions and we 
omit the proof here. ■ 

In the sequel, we expand the notation yx in Eq.Q to 
any probability vectors. Another equivalent expression of this 



lemma is that x G S^{y) if and only if x ^ yx- In this 
paper, we will switch between these two expressions from time 
to time for convenience. This simple lemma is quite useful 
because it establishes a relationship between probabilistic 
transformations and deterministic ones, while the latter have 
been well researched. 

We know that S{y) is just the convex hull of all vectors 
which may be obtained by permutating the components of y. A 
direct application of the above lemma is a similar description 
of the generalized set S'''{y), as the following theorem shows. 

Theorem 1: For all y G T/" and < A < 1, S^{y) is 
compact and convex. Furthermore, S^{y) is the convex hull 
of the vectors in the following set 

{Pyx ■ P is any n dimesional permutation}. 
Proof: Direct from Lemma^and the known structure of 

Siyx). ' ■ 

The next theorem shows that the set S^{y), as a function of 
A, is monotonic, and the intersection of all >5'^(?/), < A < 1, 
gives rise to S{y). 

Theorem 2: Suppose y G V" and < Ai < A2 < 1. Then 
S^^iy) ^ S^^{y). Furthermore, we have that 

0<A<1 

Proof: The monotonicity of S^{y) is obvious from the 
definition. So 

^(y)c fl S\y) 

0<A<1 

holds. To show 

Siy)2 n S\y), 

0<A<1 

suppose X e no<A<i ^^(y)- follows that x ^ yx, or 
equivalently, 

Eiix) > Ei{yx) ^ XEi{y) 

for all 1 < ; < n and < A < 1. When A tends to 1, we 
have Ei{x) > Ei{y) for all 1 < ^ < n. Thus x G S{y). This 
completes the proof. ■ 

From the monotonicity of S^{y) as a function of A, we 
can define the notions of limit as follows. We call S^{y) left 
continuous at A if for all nondecreasing sequences {A; : i = 
1,2,.. .}, lim^^oo Ai = A implies that 

oc 

[]S^^{y) = S\y). 

1=1 

While S^{y) is said to be right continuous at A if for all 
nonincreasing sequences {A^ : i = 1, 2, . . .}, Mmi^ao Xi — X 
implies that 

00 

\JS^^{y) = S\y). 

1=1 

Furthermore, S^{y) is continuous at A if it is both left 
continuous and right continuous. Having these notions, we are 
able to present the following theorem. 
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Theorem 3: For all y E V^, is continuous at any 

A when < A < 1. It is also right continuous at and left 
continuous at 1. 

Proof: Easy to check from the definitions. ■ 

We have examined thoroughly probabilistic transformations 
without catalysis; in the following sections, we will consider 
catalyst-assisted ones. At the end of this section, we introduce 
some lemmas that are useful for later discussion. 

Lemma 2: Given y G V"^, the function P{x — > y) is 
concave in x. 

Proof: For all x^x' G and t £ [0, 1], we have 

Eijtx + {1 - t)x') 
Ei{y) 



P(tx + (1 - t)x' ^ y) = min; 



> min; 



Eijtx) + Ei{{l - t)x') 



Ei{y) 

>tP{x^y) + {l-t)P{x' ^y). 
So P{x y) is concave in x. ■ 

The next two lemmas consider the properties of the maximal 
transformation probability P{x — > y) under the operations of 
direct summation and tensor product on its parameters x and 

y- 

Lemma 3: For all x,y £ V" and x',y' G V™', 

P{x ® x' ^ y ® y') > min{P(a; ^ y),P{x' y')}, 

where © means direct summation. In particular, P{x (S) c 
y (B c) > P{x y) for all c. 

Proof: By Vidal's formula for the probability of en- 
tanglement transformation, there exists an index I such that 
1 < Z < n + m and 



P{x ®x' ^yCBy') 



Ei{x®x') 



Ei{y®y') 

We assume that Ei{x®x') — Ei^{x) + Ei^,{x') for some < 
n and l^, < m. Notice that Ei{y ® y') < Ei^{y) + Ei^,[y') 
by definition. It follows that 



P{x®x'^y®y') > 



Ei^{x)+Ei^,{x') 
EiM + Ei,{y') 



where the second inequality follows from the following fact 

a + b b a b 

; > —, ^ — > —, for any a, b,c,d> 0. 

c + a a c a 

Thus P{x®x' -> y®y') > min{P(a; -> y),P{x' y')}. ■ 

Lemma 4: For all x,y E y" and x',y' E V"\ 

P{x x' ^ y y') > P{x y)P{x' -> y'), 

where ® means tensor product. In particular, P{x ® y ® 
c) > P{x y) for all c. 

Proof: This result is obvious from the physical meaning 
of P{x y) since the way that separately transforms x into 



y and x' into y' gives an implementation of transforming x (8) 
x' Xo y ® y'. The probability of success is the multiplication 
of the probabilities of those two transformations. That means 
P{x ® x' y (g)y') > P{x -> y)P{x' y'). 

We can, however, give a simple pure mathematical proof 
as follows. Without loss of generality, we assume that the 
components of x,x',y, and y' are nonincreasingly arranged, 
respectively. For an arbitrarily fixed integer I satisfying 1 < 
I < nm, let be the smallest index of the components of x' 
in summands of Ei{x (g) x') that have the form Xix'^, where 
1 < i < n. That is. 



Vi = min{j : Xix'^ > (x (E) 



(3) 



In case of repeated values of components of x®x', we regard 
the terms with smaller i to be included in the sum first. If the 
set in the right hand side of Eq.(|3} is empty for some i (that is, 
any term having the form xix'j, 1 < j < m, does not occur), 
then let — m+1. With these notations, we can arrange the 
summands of Ei(x (g) x') as 

n m 

El {x ® x') = '^Xi'^ x'j . 

*=1 j=ri 

By the definition of P{x' — > y'), we have X]j=r ^'j — 
P{x' ^ y') E%r, y'j for all n. Thus 

71 m 

El {x x') >p{x' -*y')Y,^^Y.y'^■ 

i=l j=ri 

Now we rearrange the summands of J27=i Sj=r y'i ^^^^ 
that 



J = l i=tj 



for some 1 < ti , . . . , tm < n+1. By the definition of P{x 
y), we have 



E^ 

i—tj 



> P{x ^ y) ^ yi 



for all tj. Thus 



Ei{x ® x') > p{x' ^ y')p{x ^y)Y.y'^Y.y^ 

j=l i=tj 

>P{x'^y')P{x^y)Ei{y®y'), 

and P{x ® x' ^ y ® y') > P{x y)P{x' — > y') from the 
arbitrariness of L ■ 



III. Catalyst-assisted Probabilistic 
Transformation 

The aim of this section is to consider the case of en- 
tanglement transformations with the aid of catalysis. First, 
we present a necessary and sufficient condition for a given 
probability vector to serve as a partial catalyst for a certain 
probabilistic transformation. 

Without loss of generality, we concentrate on catalysts with 
nonzero components, since for any probability vector c, c 
and c © have the same catalysis power in the sense that 
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in any situation, if one serves as a partial catalyst for some 
transformation, so does the other for the same transformation. 

Theorem 4: Suppose x and y are two nonincreasingly ar- 
ranged n-dimensional probability vectors, and P{x y) < 

min{a;„/y„, 1}. Let 

L = {l:l<l<n and P{x ^ y) ^ ^^}- 

Ei{y) 

Then a nonincreasingly arranged fc-dimensional probabiUty 
vector c serves as a partial catalyst for the transformation from 
X to y, that is, 

P{x ® c ^ y (g) c) > P{x — > y), 

if and only if for all ri, r2, . . . , e LU {n + 1} satisfying 
ri > . ■ ■ > Tk ^ n + 1, there exist i and j, 1 < j < i < k, 
such that 



Ci yrj Ci yr^-1 

— < — or — > — - — 

Vri — 1 Uvi 



(4) 



Here, in order to avoid a too complicated statement, we 
ignore some extreme cases of Eq.(|^; the condition © should 
be understood in the following way: whenever one of the two 
components of the disjunction in Eq.Q contains the meaning- 
less term yn+i, it is considered to be violated automatically, 
and the other component is then required. 

Proof: We will prove the theorem by showing that c 
cannot serve as a partial catalyst for transforming x into y, 
that is. 



P{x^ 



y^c) = P{x y), 



if and only if there exist ri, r2, . . . , e iU{n + l} satisfying 
ri > . . . > rfe 7^ n + 1, such that for all 1 < j < i < fc. 



Vn-l 



Vr 



(5) 



where any constraint containing the meaningless term y„+i is 
considered to be satisfied automatically. 

Notice that for all Z, 1 < ^ < nk, we can arrange the 
summands of Ei {x (g) c) such that 



Eiixi 



where 1 < r j < n+1. The case rj — n+1 means that any term 
having the form CjXi, 1 < i < n, does not occur Without loss 
of generality, we regard terms with smaller j to be included in 
the sum first in case of repeated values of components of x^c. 
This assumption guarantees that ri > . . . > r^. Furthermore, 
we exclude the possibility of ri ~ . . . — — n + 1 from 
J2j fj — I and 1 < Z < nk. 

From the definition of Ei{y (g) c) and P{x y), the 



following inequalities are easy to check: 



Eijx^c) ^ 
Ei{y®c) - k 



j = l i=rj 



> 



k n 

j=l i=rj 

k n 

j=l i=rj 



= P[x^y). 

The first equality holds if and only if Ei{y ® c) = 
Ylj^i X]r=r y*' while the second equality holds if and only 
if every is in L U {n + 1} . Consequently, we see that 
P{x c — > y ® c) = P{x y) if and only if there exist 
ri, r2, . . . , r/j G L U {n + 1} such that 

k n 

Ei{x(»c) ^'^Cj'^Xi (6) 

j=l i=rj 



and 



Ei{y ® c) Cj ^ yi 

j = l i=rj 



(7) 



for some \ <l < nk. 

In what follows, we derive the conditions presented in 
Eq.(|5} from Eqs.(|6} and 0. In fact, Eq.(0 means that 
maxi<i<fc{?/r,cj < mini<j<fe{2/ri-iCi}, or equivalently. 



Vri-l Cj 



(8) 



for all i,j = 1,2,..., A: and i > j. The special case of 
Ti = n + 1 01 rj = n + 1 can be included in Eq.(|8} by simply 
assuming that the constraints in Eq.(|8ji containing the mean- 
ingless term yn+i are automatically satisfied. Analogously, we 
can show that Eq.(|6|l is equivalent to 



(9) 



for all i > j. Notice that for all r^, G L, 

— ^ < P{x -^y)< . 

It follows that the constraints in Eq.(|9j can be derived from 
those in Eq.(|8}. That completes our proof. ■ 

Intuitively, if we decompose x (g) c and y (K) c as 



and 



X (g) C = C\X © ... © CfeX 



y ® c = ciy © . . . © Cfcy, 



respectively, then when the conditions in Eq.(|8jl are satisfied 
for some ri,r2, . . . ,rfc e L U {n + 1}, we have CiXn < 
CjXrj^x and CjX^^ < CiX^.-i for all 1 < i,j < n. So the 
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smallest k{n + 1) — X]i=i '''i components of a; ® c are exactly 
the components of the form XiCj, where 1 < j < k and 
Tj < i < n. A similar argument holds for y ® c. It follows 
that when we take I = — k + 1, then Ei{x c) — 

P{x y)Ei{y (g) c) and thus P{x ^ y^c) = P(x y). 

To our surprise, the constraints presented in Eq.Q for 
the probability vector c to serve as a partial catalyst for 
transforming x into y are almost irrelevant to x. The only 
effect of X is to determine the index set L. 

Corollary 1: Let x, y, c, and L be as in the above theorem. 

If 

P{x (8) c ^ y (g) c) > P{x y), 



then 



> max{ — : / e L\. 

Ck-i yi 



(10) 



Proof: For any I £ L, take ri = . . . = rk-i = + 1 and 
Tfe = I. Noticing that the constraints having the term y„+i are 
violated automatically, we can reduce Eq.0 to the condition 
that 

(11) 



Ck ^ yn 
Cj yi 



for some j < k. So 



Cfc-i Cj yi 
and the corollary holds from the arbitrariness of I. ■ 

Intuitively, Corollary ^ shows that the difference between 
the smallest two components of a partial catalyst cannot be too 
large. The following corollary, on the other hand, shows that 
the difference between the smallest and the largest components 
of a partial catalyst cannot be too small. Recall that a uniform 
state is a state which has equal nonzero Schmidt coefficients. 

Corollary 2: Let x, y, c, and L be as in the above theorem. 

If 

P{x (8) c ^ y (g) c) > P{x y), 



then 



Ck ^ . r yi 

ci yi-i 



I e L}. 



(12) 



In particular, any uniform state cannot serve as a partial 
catalyst for any probabilistic transformation. 

Proof: For any I £ L, let ri = ... = rk = L Then the 
conditions in Eq.0 become 



1 < JL- 

cj yi-i 



(13) 



for some i > j. Noticing that ^ < ^ for all 1 < ?' < i < /c 

Ci Cj 

and I is taken arbitrarily from L, we complete the proof. ■ 

A special and perhaps more interesting case of Theorem |3 
is when L has only one element, that is, L = {1} for some 
1 < / < n. In this case, we have the following corollary. 

Corollary 3: Let x, y, and L be as in Theorem|3] If i = {/} 
for some 1 < I < n, and 



Vn 

yi 



< 



yi 

yi^i 



(14) 



then any nonincreasingly arranged fc-dimensional probability 
vector c with 

y^<'±<JL- (15) 

yi ct y/-i 

for some t < k serves as a partial catalyst for transforming x 
into y. 

Proof: Since L = {/} contains one element, any choice 
of n, r2, . . . , rfc e _LU{n+l} satisfying ri > . . . > rk ^ n+\ 
has the form 



ri = 



Th 



1, 



Th+l 



Tk 



I 



for some Q < h < k. Take i — k and j ~ t. In what follows 
we show that under this choice of Ea. dlSt implies Eq.Q 
in Theorem |3 and thus c is a partial catalyst for transforming 
X into y. 

In fact, if t < ft,, then = I and rj = n + 1. Thus 

yr,-l _ yn ^ Cfe _ C£ 



yr. 

Similarly, if t > h then r.^ 

yr, 



yi 



Ct ( 

- I and 



yi~i 



yi ^ £fc _ 



Ct 



Thus the conditions in Eq.Q holds and that completes our 
proof. ■ 

When 2-dimensional catalysts are considered, Ea.( ll5> re- 
duces to 

^ < ^ < ^. (16) 

yi ci yi_i 

Furthermore, it is easy to check that Ea.( ll6> is indeed a nec- 
essary and sufficient condition for a 2-dimensional probability 
vetor c to be a partial catalyst. Thus Ea.(ll4> also becomes a 
necessary and sufficient one to guarantee the transformation 
from X to y has 2-dimensional partial catalysts. 

From Theorem 13 we can derive a necessary and sufficient 
condition for when a given probabilistic transformation has 
partial catalysts. 

Theorem 5: Suppose x,y £ V"^ and the components of 
X and y are nonincreasingly arranged, respectively. Then the 
probabilistic transformation from a; to y has partial catalysts 
if and only if 

P{x ^ y) < min{a:„/y„, 1}. 
Proof: The 'only if part is easy and we omit the details 
here. The proof of 'if part is as follows. 

We abbreviate P{x ^ y) to P for simplicity in this proof. 
Let l,nin and Imax be the minimal element and maximal 
element in L, respectively. Then y„ < yi^^^ since otherwise 



''max 



>P. 



yn 



which contradicts the assumption that Imax G L. Now let a be 
a real number such that yn/yimax < ct < 1 and k a positive 
integer such that a'^^^ < yi„„^/y/„i„-i- In what follows, we 
prove that 



(1, a, 



,a 
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will serve as a partial catalyst for the transformation from x 
to y, that is, 

P{x ® y® c) > P{x — !■ y). 

Here we omit the normalization of c for simplicity. 

For all ri,r2, ■ ■ ■ ,rk G L U {n + 1} satisfying ?'!>...> 
Tk ^ n + 1. There are two cases to consider 

Case I. ri E L. In this case, let i = fc and j = 1. Then 

fi ^ ^k-l ^ j/Uax < Vri ^ Vr, 

Case 2. ri = 71 + 1. In this case, denote by to, 1 < to < 
k—1, the (unique) integer such that r„i = ri + 1 but r^+i G i- 
Now let i = TO + 1 and j = m. Then 

c,; ^ y™ ^ Vn Vtj-I 
— = a > > = . 

VI, yr,„ + i Vr, 

In either case, the constraints in Eq.Q are satisfied. So from 
Theorem 13 we know that the present theorem holds. ■ 

To illustrate the utility of Theorem |5l let us see an ex- 
ample from [11] (it has been presented in Introduction). Let 
X = (0.6,0.2,0.2) and y = (0.5,0.4,0.1). We show how to 
construct a partial catalyst by the above theorem. It is easy to 
check that L = {2} and 

— = 0.25 < 0.8 = — . 

2/2 Vi 

So we can take k = 2 and any two dimensional nonnormalized 
vector c = (l,a) for 0.25 < a < 0.8 serves as a partial 
catalyst for transforming x into y. Furthermore, from the 
remark behind Corollary |3] {c = (1,Q;) : 0.25 < a < 0.8} is 
exactly the set of all two dimensional partial catalysts for this 
transformation. 

Suppose now x is just as above while y = (0.5, 0.3, 0.2). 
Then L = {2} but yg/ya = 2/3 > ya/j/i = 0.6. So any 
two-dimensional state cannot serve as a partial catalyst for the 
probabilistic transformation from x to y. We can, however, 
construct a higher dimensional partial catalyst from Theorem 
|5]as follows. First, take a real number a > 2/3/2/2 = 2/3. In 
order not to make k too large, we should take a as small as 
possible. For example, a = 0.67. Then from the constraint 
Q,fc-i ^ ^2/2/1 = 0.6 in the theorem, we have fc > 3. Thus 
the state c = (1,Q!, a^) can serve as a partial catalyst for 
transforming x into y. 

It is worth noting that the catalyst c presented in the proof 
of the above theorem can be replaced by a sequence of 2- 
dimensional vectors. To see this, notice that the only constraint 
on the dimension k of the catalyst c is a''^^ < yi,„„^/2/;,„i„-i> 
that is, for all sufficiently large fc, 

c = (1,0;, a^, . . . , a'^^^) 

is an appropriate partial catalyst which can increase the 
maximal probability of transforming x into y. In particular, 
take k — 2'"+^ — 1 for some positive integer to. From the 
simple fact that 

(1, a, a^, . . . , a^"^'"^) = (1, a) «) (1, a^) (g) . . . (g) (1, a^"), 



the effect of the catalyst in the left hand side can be imple- 
mented by the sequence of 2-dimensional catalysts listed in 
the right hand side. From this observation, a potential 'cata- 
lyst bank' need only prepare sufficiently many 2-dimensional 
catalysts in order to help probabilistic transformation. 

We state the arguments above as the following theorem. 

Theorem 6: The set V"^ constitutes a complete set of partial 
catalysts for all probabilistic entanglement transformations. 
That is, for all positive n and x,y £ V''\ if P{x (g c ^ 
y ® c) > P{x y) for some c, then there exists a sequence 
of probability vectors ci , . . . , c,„ G such that 

P{x g) ci g) . . . (K) Cm ^ y (g) ci (g) . . . g) Cm) > P{x — > y) . 

We have presented a necessary and sufficient condition 
under which a given entanglement transformation has par- 
tial catalysts. Furthermore, the proof process constructs real 
catalysts. The constructed catalysts are, however, not very 
economical in the sense that they are usually not with the 
minimal dimension among all the probability vectors which 
can serve as partial catalysts for this transformation. In what 
follows, we show how to construct economical ones. 

First, from Theorem|3 we can design an efficient algorithm 
to decide whether a probabilistic transformation has partial 
catalysts with a given dimension. 

Theorem 7: Suppose a;, y are two n-dimensional probabil- 
ity vectors and P{x y) < min{a;4/2/i, !}• Let fc be a given 
positive integer. Then the problem of whether there exists a 
fc-dimensional partial catalyst c for transforming x into y can 
be decided in polynomial time about n. 

Proof: Without loss of generality, we assume that the 
components of x, y, and c are respectively arranged non- 
increasingly. Notice that from the proof of Theorem |4] (see 
Eq.©), the necessary and sufficient condition under which c 
can increase the maximal probability of transforming x into 
y can be reexpressed as, for all ri, r2, . . . , rfc G LU {n + 1} 
satisfying ri > . . . > Vk n + 1, 

max {ynCi} > min {yr,-iCi}. 

l<i<k l<i<k 

This condition leads to the following algorithm to decide 
whether a fc-dimensional partial catalyst exists for transform- 
ing X into y: 

1. Calculate P{x — > y) and determine the set L = : 1 < 
Kn and P{x^y)^^}. 

2. For all fc positive integers ri,r2, ■ ■ ■ ,rk chosen from LU 
{n + 1}, if ri >...> rk ^ n + 1, then solve the following 
inequality about ci, . . . , c^: 

max {2/r,Ci} > min {y^.-icj. (17) 

l<i<k l<i<k 

Then there exists a fc-dimensional partial catalyst c if and only 
if the intersection of the solution areas of Eqs. in ( I17> is not 
empty when ri , r2, . . . , range over LU {n + 1} but satisfy 
ri > . . . > rk n + 1. Notice that the solution area of Ea. fTTb 
for a given sequence ri, r2, . . . , is just the union of those 
of the following fc^ inequalities: 

2/r^ ^ Vrj — l^j: 

which can be solved in polynomial time of fc. Furthermore, 
the number of choices of ri,r2, ■ ■ ■ ,rk is less than (#L + 
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1)*^, where #L denotes the number of elements in L and we have x^c^ (y (8) c)x, which impHes that x E T^{y). 



obviously, #_L < n — 1. So the algorithm presented above 
runs in 0{k){^L+ 1)*^ = 0{n^) time, which is a polynomial 
of n when k is treated as a constant. ■ 

Notice that in [16], Sun et al have presented a polynomial 
time algorithm to decide whether a given entanglement trans- 
formation has fc-dimensional catalysts, that is, fc-dimensional 
partial catalysts which can increase the maximal transforming 
probability to 1. So a little modification of Sun's algorithm 
can also be used to determine whether or not a fc-dimensional 
partial catalyst exists. What we would like to point out is 
that the complexity of Sun's algorithm is 0{'n?^) while our 
algorithm presented in Theorem is 0{{^L)^). Although 
they are both exponential of fc and in the worst case #L = 
n — 2, in practice our algorithm is more efficient since is 
generally much less than n. 

Theorem Q and Theorem |5] together give a method for find- 
ing out the most economical catalysts for a given entanglement 
transformation as follows. First, we use Theorem |5l to decide 
whether partial catalysts exist for this transformation and an 
upper bound m on the dimensions of the most economical 
ones can also be derived. Second, for k = 2,3, ...,m we 
use the algorithm presented in Theorem to decide whether 
there exist fc-dimensional partial catalyst. Moreover, from the 
algorithm presented in Theorem^ the most economical partial 
catalysts can be constructed explicitly. 

IV. Structure of Catalyst-assisted Probabilistic 
Transformation 

In this section, we investigate the mathematical structure 
of catalyst-assisted probabilistic entanglement transformation. 
Given a probability vector y £ V"^ and < A < 1, denote 
by T^{y) the set of probability vectors which, with the aid 
of some catalyst, can be transformed into y with a probability 
not less than A, that is 



T^{y) = {x e V"" : P{x < 



y (Xi c) > A for some 



probability vector c}. 

The special case of A = 1 corresponds to T{y), which is 
just the set of all probability vectors that can be transformed 
deterministically into y by catalyst-assisted LOCC (for the 
definition, we refer to [12] or [17]). It is easy to check that 
the set T^{y), as a function of A, is monotonic. 

Recall that S^{y) is just equal to S{y\) from Lemma 
[2 One may wonder if T^{y) — T{y\), or if there exists 
a simple connection between catalyst-assisted probabilistic 
transformation and catalyst-assisted deterministic one. If so, 
then all the known properties of T{y\) can be used to give 
simple proofs of those of T^{y). In fact, we can prove the 
following. 

Lemma 5: For x,y £ V", if x G T{y\), then x £ T^{y). 
That is, T{y^) C T^{y). 

Proof: Suppose x £ T{y\). By definition, there exists a 
probability vector c such that x ® c < yx ® c. Noticing that 

yx®c = {l-\E2{y))c®\y2c® ...(BXynC^ {y®c)x, 



But unfortunately, T{yx) 7^ T^{y). Moreover, we can show 
informally that for all z £ V", T{z) ^ T^{y) as follows. 
From Corollary 5 in [12], the boundary points of T{z) is the 
set {x £ T{z) : x\ = z[ or a;^ — z^}, while from Theorem 
|9]in this section, the boundary points of T^{y) is the set {x £ 
T^{y) : x^ = Ay^}. It is obvious that these two sets cannot 
be equal. So T{z) ^ T^iy). Since T"{y) = and T^{y) = 
T{y), we assume from now on that < A < 1. The next two 
lemmas are useful for latter discussion. 

Lemma 6: Let x,y £ V" and < A < 1. Suppose x and y 
can be decomposed into two parts respectively as 

X = x' ® Xz and y — y' ® z 

for some z such that z| < y[. Let /i denote the sum of the 
components of z. Then x £ S^{y) if and only if 



where 



1 — /iA 1 — A* 

^, ^ A - ^A 



(18) 



Proof: Without loss of generality, we can assume that 

v' = {yii,yi2,- ■ ■ 

and the components of y' have been arranged nonincreasingly. 
If X e S'^'ly), then x < yx from Lemma ^ Furthermore, we 
have yx — y" ® Az and then 



X ^ y 



where 



1 — /iA 1 — /iA 



y" = {\-\E2{y),\y^„---,Xv.J- 



A simple calculation shows that 

y" , y' 



-)> 



1 — /xA 1 — 

so the 'only if part of the lemma holds. 

The 'if part can be proved by simply retracing the argu- 
ments above. ■ 

Lemma 7: Let x,y,X,ii, and A' be as in Lemma |6l Then 
X £ T^{y) if and only if 

1 — /iA 1 ^ 

Proof: Notice that x £ T^iy) if and only if there exists 
a probability vector c such that x^c £ S^{y(E)c). The current 
lemma is just a simple application of Lemma |6l ■ 
The following theorem shows that some properties of T{y) 
such as convexity and containing corresponding no-catalysis 
case as a subset are also shared by T^(j/). 
Theorem 8: For all ?/ e and < A < 1, 

1) T'^{y) is convex; 

2) S\y) C T\y)- 

3) suppose x £ T^{y) and the components of x and y 
are nonincreasingly arranged respectively. If x 7^ yx, then 
^d/yd > A, where d is the maximal index of the components 
such that Xd/yd^ A. 
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Proof: 1) Suppose x,x' E T^{y). By definition, there 
exist probability vectors c and c' such that P{x®c y®c) > 
A and P{x' ®c' ^ y®c')>\. For alH, < i < 1, we have 

P{(tx + (1 - t)x')®c^ y®c) 

> tP{x®c^ y(g)c) + (1 -t)P{x' y®c) 

> tP{x (g)c-^y(g)c) + (l- (8) c' ^ y (g) c') > A 

from Lemma 12 and Lemma 0] where c = c'. 

2) It is obvious from Lemma |4] 

3) Suppose a; G T^{y). Decompose x as a; = x' ®x" where 
x' ^ {xi,X2, ■ ■ ■ ,Xd) and x" ^ {xd+i, ■ ■ ■ ,Xn)- 

The vector y is decomposed analogously. From the definition 
of d, we have x" = Ay". Applying Lemma we have 



(19) 



1 — /iA ' 1 — /i ' 

where /i is the sum of the components of y" and A' is 
defined by Ea. (ll8> . We can then deduce that Xd/yd > A since 
otherwise (notice that Xd/yd 7^ A by assumption) 



p{- 



^ (1 - ^J)Xd ^ y 

(1 - ^X)yd 

for any probability vector c, where k is the dimension of c. 
This contradicts Eq.(ll9>. ■ 

Notice that in 3), no constraints on the largest components of 
X and y are needed for x G T^{y), in contrast to the necessary 
condition that xi < j/i of a; G T{y). This is due to the 
asymmetry of roles of the largest and the smallest components 
in determining the maximal transforming probability. 

Lemma 8: Let y e and < A < L For all x e S^{y), 
if P{x —>?/)< x^/y^, then x is in the interior of T^{y). 

Proof: Suppose x e S^{y) and P{x —^y)< x^^/y]^. If 
P{x y) = 1, then x is in the interior of S^{y) for A < 1. 
Thus a; is also an interior point of (y) since S^{y) C T^[y). 
When P{x y) < \, by Theorem |5] there exists a partial 
catalyst c such that P{x®c y®c) > P{x ^ y) > A. Thus 
X is an interior point of T^{y). ■ 

Theorem 9: Let y e < A < 1 and x G T^(y). Then 
X is on the boundary of T^{y) if and only if x^/y^ — A. 

Proof: Notice first that x S T^{y) implies x^/y^ > A 
from Theorem|8] Suppose now that x^^/y}^ = A. For all e > 0, 
consider the probability vector 

X ~ (xi, . . . , xl^_2, Xj^_-j^ + e, xl;^ — e). 



(1 - fl)x'^Ck 



It is easy to check that 



P{x' (8) c ^ y (g) c) < 



Vn yn 



for any probability vector c, where k is the dimension of c. 
Thus x' ^ T^iy). It follows that a; is a boundary point of 
T\y). 



Conversely, suppose x\-^/y*^ > A. By the definition of x £ 
T^{y), there exists a probability vector c such that P{x®c 
y(E)c) > A. If the inequality holds strictly, then x is of course in 
the interior of T^{y) by the continuum of P{x®c y®c) on 
X. So we need only consider the case of P{x®c y®c) = A. 
In this case we have 



P{x ®c^y®c) = \< 



where k is the dimension of c. Thus x ®c is an interior point 
of (y (g) c) from Lemma |8] On the other hand, since the 
function x i-^ a: (g c is continuous, it follows that x is in the 



interior of the set {a 



c G T^{y ® c)}, which is just a 



subset of r (y). So x is in the interior of T [y). ■ 

Recall that the boundary point set of T{y) is {x € T{y) : 
— y\ or x\^ = y/j}. Once again, the asymmetry of roles 
of the largest and the smallest components in determining the 
maximal transforming probability makes the boundary point 
set of T^{y) different from that of T(y). 

The next theorem tells us when catalysis is helpful for 
probabilistic transformation with destination state y by giving 
a necessary and sufficient condition for T^{y) = S^{y). 

Theorem 10: Let y £ V" be a nonincreasingly arranged 
probability vector and < A < 1. Then T^{y) = S^{y) if 
and only if y2 = yn- 

Proof: If y2 = yn, then for any nonincreasingly arranged 
X E and any integer / satisfying 1 < I < n, 



Ei[x) 



En 
i=l ■ 



> 



{n ~ I + l)Xn 



Xn 
Vn 



Ei{y) {n~l + l)yn (n - Z + l)y„ 

Thus P{x y) = min{x„/y„, 1} and for any vector c, P{x® 
y®c) = P{x y). It follows that T^(y) = S^{y). 
Conversely, suppose y2 > yn- Let m > 1 be the maximal 
index of the components of y which are not equal to y„, that 
is, yrn > IJrn+i = -■ - = yn- Then we have En,{y) > (n - 
m + l)y„- Let /i be a real number such that 

A£^m(y) 



X < fi < min{ 



(n 



,1}, 



and define a probability vector 

1 - A 



(2/1 



1 



£',„(y), . . . ,ym_i + 



1-A 



Ai^m(y) - M£^m+i(y),Mym+i, . 



1 



It is a little tedious but very easy to check that the components 
of X have been nonincreasingly arranged and P{x ^ y) = 
A < fj, = Xn/yn- By Theorem |9j x is an interior point of 
T^{y). That completes our proof that S^{y) 7^ T^{y) since x 
is obviously on the boundary of S^{y)- ■ 

Compared with the conditions of T{y) = S{y) presented in 
[12], our condition in Theorem 1101 is more simple and even 
a little surprising since it depends only on whether or not 
2/2 — Vn (totally irrelevant to the value of A). 

In what follows, we consider the interesting question of 
whether or not there exists a bound on the dimension of partial 
catalysts that we should consider to help transforming some 
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vector X into a given probability vector y. This is in fact a 
generalization of the problem considered for T{y) in [12]. We 
will give a negative answer to this question by showing that 
in general T^{y) 7^ T^{y) for all positive k, where 

Tkiy) = {xeV" : P(x®c y(»c) > A for some c e V''} 

is defined to be the set of probability vectors which, with the 
aid of some fc-dimensional catalyst vector, can be transformed 
into y with the maximal probability not less than A. 

Lemma 9: For all y e V"" and any positive integer k, [y) 
is a closed set. 

Proof: Suppose a;^, x^, ... is an arbitrary vector sequence 
in T^{y) that converges to x. By the definition of T^{y), there 
exists a catalyst sequence c^, c^, . . . in V'^ such that P(x' (g) 
c* y®d) > A for i — 1,2, . . .. Notice that V'' is a compact 
set in R'^. There exists a convergent subsequence d^, . . . 
of c^,c^,... that converges to, say, c E V''. Then we can 
deduce that P{x ^ c ^ y ^ c) > A and so a; G T^{y) from 
the fact that the function P{x ® c y ® c) is continuous on 
the parameters x and c. ■ 

Notice that in [12], a similar lemma about Tk{y) was 
presented but the proof there was a little complex. The proof 
technique of the above lemma can be used to give a simpler 
one. 

Theorem 11: Let y E be a nonincreasingly arranged 
probability vector, ?/2 > Un, and < A < 1. Then T^{y) ^ 
T^{y) for all positive k. 

Proof: We will complete the proof by showing that under 
the assumptions in this theorem, T^{y) is not a closed set. 
Then from Lemma|9] we have T^{y) 7^ T^{y) for all positive 
k. 

Let e — {1/n, . . . ,l/n)he the uniform vector in V". Notice 
that e E T^{y) and from Theorem|9l e is also an interior point 
of T^{y). Denote by m the maximal index of the components 
of y which are not equal to j/„, that is, 

Vm > Vrn+l =...=?/„. (20) 

Since ?/2 > y„, we have 1 < m < n. Let fj, E (0, A) and 
define a probability vector 



1 — /i 

a; = (2/1 H -E,n{y), ym-i 



1 



1 



E,n{y), 



It is easy to check that the components of x are arranged 
nonincreasingly and P{x y) — a;„/y„ — jj, < X. So x ^ 
T^iy). If T^{y) is closed, then the closed set {tx + (1 - t)e : 
< t < 1} must intersect T^{y) at a boundary point of T^{y), 
say, x' = t'x + (1 — t')e. Thus x'^jyn ~ A from Theorem |9l 
On the other hand, from Ea. (l20> and noticing that < < 1, 
we have 



Vn 



\-t' 

nyn 



l-t' 

ny„i 



x^ 

Vn 



and so 



■'m+l 



y-m Vm+l 



Xn 
Vn 



Notice that when y2 = yn, we have T^{y) = S^{y) from 
Theorem 1 101 and so T^{y) = T^{y) for any positive integer 
k since 

T\y) = S\y) C T^{y) C T\y). 

In order to describe T^{y) more deeply, we examine its 
extreme points in what follows. 

The following lemma describes what kind of perturbations 
will not remove a point from T^{y), even when the point 
is on the boundary. This lemma can also be regarded as a 
generalization of Corollary 5 in [12]. 

Lemma 10: Suppose < A < 1, and x,y E be two 
nonincreasingly arranged n-dimensional probability vectors 
satisfying x E T^{y) but x ^ y\. Let d be the maximal 
index of the components such that Xd > Xyd- Then any 
sufficiently small perturbation will not remove x from T^{y), 
provided that the perturbation does not affect the components 
Xd+i, ■ ■ ■ ,Xn- (Notice that if d = n, then x is an interior 
point of T^{y) and the result here is just a simple property of 
interior points.) 

Proof: To begin with, let us decompose x as x — x' ® x" 
where 

x ^ (xx,X2,. ■ ■ ,Xd) and x" = (x^+i, . . . , x„). 

Similarly, y can be decomposed as y = ?/' y" . Applying 
Lemma we have 



xET\y) ^ xET^'iy), 



where 



X' y 
X = r and y = 



(21) 



1-AiA 



and /i is the sum of the components of y" and A' is defined 
by Eq.(ll8>. Since Xd > Xyd, we have 

Xd ^ w yd w~ 

Xd = r > A = A yd, 

1 — flA 1 ^ M 

and so x is an interior point of T^{y) from Theorem |9] 
Notice that any perturbation to x which does not affect the 
components Xd+i, ■ ■ ■ , Xn is also a perturbation to x and vice 
versa. Therefore the lemma follows. ■ 



Using this lemma, we can easily 
points of T^{y) for a given y E V 
that T^{y) and S^{y) share the same 
they satisfy very different properties 
in general. Notice that T^{y) is not 
proof of Theorem so we cannot 
T^{y) only by its extreme points. 

Theorem 12: For all y E V", the 
T^{y) is the same as that of S^{y). 



find out all the extreme 
"■. It is rather surprising 
extreme points, although 
and also S^{y) g T^{y) 
closed in general by the 
determine the whole set 

set of extreme points of 
That is, it is just the set 



{Pyx '■ P is any n dimesional permutation}. 
Proof: First, we prove that any extreme point of S^{y) is 
an extreme point of T^{y). To prove this, we need only show 
that y\ is an extreme point of T^{y). Suppose the components 
of y have been arranged nonincreasingly and 



which contradicts 3) of Theorem |8l 



yx = ty' + (1 - t)y" 



(22) 
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is a convex combination of y' and y" in T^{y), where < 
t < 1. From 3) of Theorem |8] we have 

y'n > y'n > Ay„ 

and 

Vn > y'n^ > Ayn 

(notice that y' and y" are not necessarily nonincreasingly 
arranged). But from Ea.( l22> we have ty'„ + (1 ~ t)y'/^ — Ay„. 
It follows that 

y'n = y'n = Ay„, (23) 

and furthermore, y'^ and y'^ are the smallest components of y' 
and y", respectively. Now using 3) of Theorem |8] again, we 
can deduce that 

y'n-l y'n-1 A2/„_l 

from Eas. (l23> and ( I22> . And again, y'^^i and are the sec- 
ond smallest components of y' and y", respectively. Repeating 
the above arguments, we finally come to the conclusion that 

y' = y" = yx- 

Thus y\ is an extreme point of T^{y) as claimed. 

Now suppose X is an extreme point of T^{y) and x is 
nonincreasingly arranged. If it is not an extreme point of 
S^{y), then x ^ y\ and so there exists a positive integer 
d, 1 < d < n, such that Xd > Xyd- By Lemma [Tol we can 
find a sufficiently small but positive e such that x' e T^{y) 
and x" G r^(?/), where 

x' ^ {xi, . . . , Xd-2, Xd-i - e,Xd + e, Xd+i, ■ ■ ■ , Xn) 

and 

x" = {xi, . . . , Xd~2, Xd-1 + e,Xd- e, Xd+i, ■ ■ ■ , Xn)- 

Obviously x = (x' + x")/2, which contradicts our assumption 
that X is an extreme point of T^{y). That completes our 
proof. ■ 
What we would like to point out here is that a similar 
argument on T{y) instead of T^{y) in the above theorem can 
lead to a solution to the open problem Nielsen proposed in his 
lecture notes on the theory of majorization and its applications 
in quantum information theory [14]. More specifically, T{y) 
has a discrete set, but not a continuum as Nielsen conjectured, 
of extreme points (in fact, the extreme points of T{y) are just 
{Py : where P is any n-dimensional permutation}). 

V. Conclusion and Open problems 

In this paper, we investigate carefully the power of catalysis 
in probabilistic entanglement transformations by LOCC. We 
give a necessary and sufficient condition for when some 
appropriately chosen catalyst can be helpful in increasing the 
maximal probability for a given probabilistic transformation. 
An efficient algorithm is presented to decide whether partial 
catalysts with a given dimension exist for a certain probabilis- 
tic transformation, which leads to a method for constructing 
the most economical partial catalysts with minimal dimension. 
We also study the set of states that can be transformed by 
catalyst-assisted LOCC into a given state with the maximal 



probability not less than a given positive number. We prove 
that this set shares many properties with the well known set 
consisting of all vectors being trumped by the given state. 
More mathematical structure of catalyst-assisted probabiUstic 
transformation is also considered. 

The emphasis of this paper is to determine when the 
maximal transforming probability can be increased in the 
presence of partial catalysts and how to construct appropriate 
ones. The amount of the probability increased is, however, 
not considered. So an open problem and also an important 
direction for further study is to determine the maximum of the 
transforming probability that can be reached with the aid of 
partial catalysts. This is also a generalization of deterministic 
catalysis since the problem of the existence of catalysts for 
deterministic transformations is equivalent to the problem 
of the existence of partial catalysts which can increase the 
transforming probability to 1 for probabilistic transformations. 

At the end of Section IV, we have determined all the extreme 
points of T^{y). However, we still know little about the 
geometric structure of T^{y). The main reason is that T^{y) 
is not closed in general and so it is not the convex hull of 
its extreme points. How to determine the accumulation points 
outside T^{y) is really a hard problem and remains open. We 
believe that T (y), the closure of T^{y), has a continuum of 
extreme points, just as Nielsen conjectured. 
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